Local Inhomogeneity Effects on Nucleation Process 
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Quantum nucleation processes in the presence of local 
moderate inhomogeneities are studied theoretically at high 
biases. The quantum nucleation rate Gamma is calculated for 
one- dimensional systems in a form V = Ae'^ — B /K) by using 
the 'bounce' method. The bias-dependence of the exponent 
B is shown to be changed by inhomogeneities. This change is 
explained by the reduction of the effective spatial dimension 
of the system. By studying the system-size dependence of the 
prefactor A, the condition for the appearance of inhomogene- 
ity effects is evaluated. Nucleation rates in thermal activation 
regimes are also calculated, and compared with quantum tun- 
neling regimes. For higher-dimensional systems, it is shown 
that the local approximation of inhomogeneity does not hold, 
and that spatial profiles of inhomogeneity become important. 

FAGS numbers: 64.60. Qb, 03.65.Sq, 73.40. Gk 



I. INTRODUCTION 

Nucleation process is one of the most universal phe- 
nomena found in various areas in physics from cosmol- 
ogy to condensed matter physics providing the mecha- 
nism for the onset of first-order transitions Besides 
classical nucleation due to thermal fluctuations, nucle- 
ation due to quantum tunneling has attracted interest 
for many years as macroscopic quantum phenomena. Al- 
though quantum nucleation has been studied theoreti- 
cally in the pioneering issue 1^, it is just recently that 
observation of quantum nucleation has become possible 
owing to the progress in experimental techniques. At 
present, quantum nucleation has been observed experi- 
mentally in low-temperature condensed-matter systems: 
^He-'^He liquid solutions [pP], cavitation in ''He liquid 
1^, nucleation of ''He solids^. Nucleation of magnetic 
domain in thin films has also been discussed theoreti- 
cally [0-^ . Quantum creation of a kink-antikink pair can 
be regarded as one-dimensional nucleation, and has been 
studied both experimentally and theoretically in disloca- 
tion motion in solids JPo 11 1 , and sliding of charge density 
waves (CDWs) |17]], though in the latter system in- 
terpretations of the experimental results are not settled 
I believe that long Josephson junctions 21-2^] 



are also suitable to observe quantum nucleation of a soli- 
ton pair. 

In many theoretical issues, nucleation processes have 
been restricted to the situation that a stable phase is 
formed in the otherwise homogeneous background of an 
unstable phase. In several literatures, the nucleation 
rates in the presence of strong inhomogeneities have also 
been studied by using a single- variable model |26|. In 
such studies, direct evaluation of the impurity strength 
from the nucleation rate is difficult, because the parame- 
ters of the single- variable model cannot be related clearly 
to the impurity strength. In this paper, we study mod- 
erate inhomogeneities, which can be controlled by an ex- 
ternal parameter from the nearly-homogeneous region to 
the strongly-inhomogeneous region. This crossover re- 
gion has not been studied in detail, and provides useful 
information of the impurity strength, because the nucle- 
ation rates are very sensitive to the impurity strength in 
this region. 

In this paper, we study how the nucleation process 
changes in the presence of local moderate inhomogene- 
ity. As a starting point, a highly-biased region is inves- 
tigated where the potential can be expressed by a cu- 
bic polynomial. The nucleation rate is calculated based 
on the 'bounce' method p^^-|29[| for quantum tunneling 
regimes, and on Kramers' law for thermal activation 
regimes [|l|,Q . We concentrate on nucleation processes in 
one-dimensional systems which correspond to the kink- 
antikink nucleation. We show that extension to higher 
dimensions is not easy since the nucleation rates are af- 
fected by details of local inhomogeneity profiles such as 
an impurity size. 

The nucleation rate is expressed in the form F = 
Aexp(— _B/?i), where A and B are called as a prefactor 
and an exponent, respectively. The exponent B depends 
on the bias / as 



B^Uc- fV 



(1.1) 
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where fc is a classical threshold bias at which the po- 
tential barrier disappears. It should be stressed that the 
value of 7 is modified by inhomogeneities. The mod- 
ification of 7 has first been discussed in the quantum 
sHding of CDWs by Yumoto et.al (3|-|3|. They have 
calculated the nucleation rate by reducing the model to 
a single- variable problem with use of path integrals, and 
have reported that the value of 7 increases as the im- 
purity strength is enlarged. This method, however, is 
so complicated that it is difficult to discuss the origin of 
the change of 7 clearly, and also to examine the valid- 
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ity of the approximation adopted there. In this paper, 
we study inhomogeneity effects without any reduction to 
a single- variable problem. It is claimed that our results 
are explained by the 'dimensionality' of nucleation, which 
has been pointed out by the author in ref. |Q . 

We consider one isolated impurity in a system with the 
size L. The results obtained in this paper are also ap- 
plicable to systems with the dilute impurities by taking 
the average impurity distance as L. In the limit L ^ oo, 
nucleation occurs dominantly in homogeneous regions of 
the samples, and no impurity effects appear. The inho- 
mogeneity effects appear only when L is below a crossover 
value Lcr. Within the bounce method, Lcr is shown to be 
obtained by calculating the prefactor A. Here, we should 
note that the validity of the bounce method is not guar- 
anteed generally for the many-body problems. We also 
discuss the validity of the bounce method when it is ap- 
plied to the nucleation problem in the presence of the 
impurities. 

This paper is organized as follows. The model Hamil- 
tonian is given in Sec. |ll[ The quantum tunneling rate is 
formulated and calculated for one-dimensional systems 
in Sec 



III 



IV 



We study the thermal regime in Sec 
briefly, and the extension to higher-dimensional systems 
and justification of the bounce method are discussed in 
Sec. |V|. Finally, results are summarized in Sec. VI. 



Another example is the (j)'^-mode\ with a bias 



(2.4) 



The type in ( |2.3| ) has been used for long Josephson junc- 
tions, dislocation in solids, and CDW systems, while the 
type in ( p.4[ ) describes an effective model for nucleation 
of a stable phase from an unstable background. 

In this paper, the characteristic action of the systems 
is assumed to be large as compared with the Planck con- 
stant. In such situations, the quantum nucleation ap- 
pears only when the bias / is controlled near the classical 
threshold fc at which the metastable state of the poten- 
tial Vo{(j)) disappears. Note that /c,o = 1 for the tilted 
cosine potential and /c,o = 2\/2/3 for the i/i'^-model in 
the abs ence o f im purities. Near / = fc^, the potentials 
in ( ^.3|) and (2.4) form can be approximated commonly 
by quadratic-plus-cubic polynomials around a metastable 
point (p = (po as 
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(0-0o)' 



(2.5) 



II. MODEL 



Here, a = 0(1) and b = 0(1). To be explicit, a = 2^/^ 
6=1 for the tilted cosine potential (2.3), while a — 2"^/^, 
b = 2^/2 for the ^^-model (O. 



In this paper, nucleation rates are calculated based on 
the equation of motion 

dV 

- V^^ -f — = 0. (2.1) 

Here, (j){x,t) is a (d -I- l)-dimensional field, and x and 
t are scaled by the characteristic length and time, re- 
spectively. Dissipation is assumed to be weak enough, 
but not extremely weak so that the system is in a ther- 
mal equilibrium in a metastable well. We concentrate on 
one-dimensional systems described by the equation 

dV 

^tt - cl)xx + -g^ = 0, (2.2) 

and the higher-dimensional systems are discussed only in 
Sec. ^ The potential energy is assumed to consist of 
two parts as F(0, x) = Vb((/i) -f Vinip{(f>,x), where Vb((/>) 
is a homogeneous part and Vi,iip(0, a;) describes an inho- 
mogeneity. 

As for the homogeneous part Vb(0), we assume: (1) 
the potential has at least one metastable point (po; (2) 
the energy barrier in Vb(0) is controlled by an external 
parameter like an external field. One typical example for 
Vo{(f>) is a tilted cosine potential 

Vb(</') = (l-cos0)-/0. (2.3) 




FIG. 1. String motion in the potential Vo (<?!>). The impurity 
effect is described by a local force for this string at x = 0. 

The equation of motion ( ^.2[ ) describes dynamics 
of a string in the potential V{(f>) shown in Fig. |^. 
When the string is initially located in a metastable 
well, it stays there for a while, and then begins to 
move towards lower energy states. If the system size 
is large enough, the transition from metastable states 
to moving states occurs through local deformation of 
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the string accompanied with the creation of one kink- 
antikink pair. This transition can be regarded as 'nucle- 
ation' when the inetastable(stable) state is related to the 
metastable(stable) phase. 

The local inhomogeneity is introduced by Vinip((/>) as 



Vinip{(l>,x) = 5(x)h[(j)), 



(2.6) 



where 5{x) is a delta function, and h{(j)) is a function de- 
fined around the metastable point. This type of impurity 
potentials has been studied theoretically in CDW sys- 
tems pl|-|33| and long Josephson junctions ||35| , ^ . The 
potential form can be obtained not only for im- 



purities but also for inhomogeneities of external forces, 
modulation of barrier heights, and edge effects. (Details 
are summarized in Appendix |^.) At high biases, since 
the impurity effect is determined by a small modulation 
of (j) around the metastable point (pQ, the inhomogeneity 
potential can be approximated as 



Mmp(</', x) - -e5{x){(l) - 4)o). 



(2.7) 



Here, e describes the impurity strength. This local po- 
tential can be regarded as a local force on a string at 
X = as shown in Fig. |[ For e > (e < 0), the string is 
attracted towards the positive (negative) direction of (f>. 
As a result, the nucleation is enhanced around x ~ for 
e > 0, and suppressed for e < 0. 

The appearance of the impurity effect is explained in- 
tuitively as follows. The nucleation rate F is expected to 
consist of two parts as F = Fbuik + Timp- Here, Fbuik is 
the nucleation rate in the homogeneous region far from 
an impurity, and Fjmp is a contribution near the impu- 
rity. Using the system size L, the bulk part is estimated 
as Fbuik ~ ^qL, where Fq is the nucleation rate per 
unit length in the absence of the impurities. The im- 
purity part Fimp is independent of L, and is estimated 
Fo exp(Af7/?ici;o), where AU is the energy 



imp 



as F 

gain at the impurity site, and ojo is an attempt frequency 
around the metastable state. If the system size L is large 
enough, the bulk part Fbuik dominates the impurity part 
Fimp, and no impurity effects appear. However, by con- 
trolling AU, Finip can dominate Fbuik- The inequality 
^ imp S> Fbuik leads 



AU > TiujologL. 



(2.8) 



Since the condition (2.8) shows of a weak system-size de- 
pendence of the condition (2.8) through the logarithmic 
function of L, it is expected t o b e feasible in experimen- 
tal systems. The inequali ty (^.8[ ), which can be derived 
more accurately in Sec. [H D , will be assumed in the 
next section. 



In this paper, we further assume e > in (2.7), because 



Fiinp never dominates the bulk part Fbuik in case of e < 0. 



III. QUANTUM TUNNELING REGIME 



A. Formulation and Scaling Analysis 

The quantum nucleation rate in the presence of the 
impurity is formulated by a semiclassical method as fol- 
lows. The partition function of the system is expressed 
in terms of path integrals as 



Z = J I?(/)(x, t) exp 



5[0(x,t)] 



rL/2 i^hp 

5[0(x,t)] = / Ax dT 

J-L/2 Jo 

+ Vo{cf>)+Vi^p{x,cl>)] 



h 

2^^ + 2^? 



(3.1) 



(3.2) 



where ?i is a dimensionless Planck constant normalized 
by the characteristic frequency, length, and energy of the 
system, t = it is an imaginary time, and f3 — \/k^T is 
an i nve rse te mpe rature. The potential forms are given 
by ( 2^ ) and ( 2/7 ) . By rescaling the variables as 



c' (ay/ fcfl - 



-1/2 
-1/2 



a\//c,o - / 



+ 00, 



(3.3) 
(3.4) 
(3.5) 



the action is rewritten as 



S[ip{x,T)] = 

^1 



fe2 



dx' 

L/2 Jo 

1 



hp 



dr' 



(3.6) 



Here, L = La^/^fcO - ff'^ and /3 = (3a^/^{U,o - ff'^ 
is a scaled length and inverse temperature, respectively. 
The impurity effect is described only through an effective 
impurity strength e defined by 



(3.7) 



This effective impurity strength depends on both the bare 
impurity strength e and external bias /. Even if e is 
fixed, we can control the strength of the impurity effect 
by changing the external current; the impurity effect is 
enhanced by controlling the external current near the 
threshold current /c. 



The partition function given in (3.1) and (3.2) is eval- 
uated by integrating the partition function around sta- 
tionary solutions up to second-order fluctuations. The 
stationary solutions are determined by 



5S 
dip 



Vx'x' +(p-^- eS{x') = 0. (3.8) 



There exist two types of solutions for this equation. One 
is a stable solution ipo{x') independent of r', and the 
other is a 'bounce' solution ipb{x',t'). Since the system 
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has a metastable state, the free energy F = —f3~^ ^ogZ 
gains an imaginary part. Then, the decay rates are 
given by r = 2IniF. This formula, derived through the 
'bounce' method, was first appHed to the thermal activa- 
tion regime |p9|,[40|, and extended to the quantum tunnel- 
ing regime |2|2j|29(] . The nucleation rate is now written 
in the form 

r = Aexp{-B/h), (3.9) 
and the exponent B and prefactor A are given as 



B = S[lpb(x, t')] - S[(po{x')], 



A 



pn l\ \ ,(0) 



-1/2 



(3.10) 
(3.11) 



Here, a(T) is a temperature-dependent factor, and takes 

1 for T < Tq, and Tq/T for T > Tq |i|, and Af^s 

(a['^''s) are the frequencies of eigenmodes around the 
bounce(stable) solution determined by the equation 



\(B,0) ,(B,0). , /N 



Note that the exponent of (/c^o — /) depends crucially 
on the spatial dimension d. The factor 3/2 comes from 
the bias dependence of the barrier height AU oc (/c.o ~ 
/)3/2, while the factor -{d+ l)/4 comes from the fact 
that the typical size of the bounce solution in the spa- 
tial(temporal) coordinates is given by (/c o — f)^^^'^ for 
each of original coordinates (x, r). (See (|3.3| ) and (3.4).) 

The impurity effect can be explained by the change in 
the dimension: the strong impurity changes the spatial 
dimensionality d to zero. The classical threshold bias 
is also modified from f^ Q to fc- As a result, the bias- 
dependence of B becomes 



S«(/c-/)*---(/c-/) 



5/4 



(3.15) 



(3.12) 



in the presence of strong impurities. In this region, the 
nucleation process is well described by a local deforma- 
tion of the field, and is reduced to a one- variable problem. 

This viewpoint is valid in the one-dimensional case as 
shown in the subsequent sections. In the d> 2 case, how- 
ever, the viewpoint of the dimensional crossover in this 
paper cannot be discussed, since the impuri ty p otential 
cannot be described by the delta function in ( ^.7[ ) , as will 
be discussed in Sec. ^ 



All the eigenvalues around the stable solution are positive 
(Xf^ > 1). The lowest eigenvalue A^^' is negative due to 
the metastability of the bounce solution ips ■ The second 
lowest eigenvalue Aj is always zero due to the trans- 
lational symmetry of tpB in the t direction. This mode 
must be treated by the Fadeev-Popov method, and the 
expression of the prefactor is modified as [pSlpql 



C. Exponent 

For one-dimensional systems, the exponent B is calcu- 
lated in a form B = aB{e). The factor a = a?{fc — f)/b^ 
describes the bias-dependence in the absence of the im- 
purities, while the scaled exponent B{£) is calculated as 



(B) 



B 

2TTh 



hp 



(3.13) 



where Tq denotes the center position of the bounce. For 

( B] 

the weak impurity, the third eigenvalue Ag also ap- 
proaches zero due to the translational symmetr y of i pB 
in the cc-direction. We study this mode in Sec. Ill D in 
detail. 



B{e) = S[Mx',T')]-S[Mx% 



S[^{x',t')] 



L/2 j-hp 

dx' / dr' 

-L/2 ^0 



2^. 



1 2 

2^ 



^ip"^ - e5{x')ip 



(3.16) 



2^^ 



(3.17) 



The impurity effect appears only through B[e). In the 
homogeneous case, B{e) is obtained as 



B. Dimensional Crossover 

The bias dependence of the exponent B can be dis- 
cussed from the viewpoint of the 'dimensional crossover'. 
This viewpoint has been discussed in several physical sys- 
tems such as CDWs [Q, and long Josephson junctions 
p7t in the context of the system-size dependences. The 
application to impurity problems has been discussed re- 
cently in ref. p^ . Here, we briefly summarize this view- 
point. 

For the weak impurity e <C 1, the bias-dependence of 
the exponent B is obtained as 



5(0) = S[ipB{x\ r'; e = 0)] 31.00, 



(3.18) 



and the exponent behaves as B — B{0)a oc {fc — fY- 
As the effective impurity strength e increases, B{£) is 
suppressed, and the tunneling rate is enhanced. At a 
critical value id — 4/\/3 ~ 2.31, the potential barrier 
disappears, and the exponent B is reduced to zero. Note 
that before the exponent B becomes zero, the bounce 
method becomes invalid because it is justified only in the 
semiclassical condition B/h ^ 1. However, the region 
where the bounce method is broken down is narrow since 
the normalized Planck constant fi is assumed to be small 
enough. 



B^{Ufl~fY' 



(3.14) 
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W e calculate the function B{e) by solving the equation 
( |3.8|) numerically. The Newton method is used by divid- 
ing the x'-t' {x' > 0, r' > 0) space into 64 x 64 cells. 
The result is shown in Fig. ^ by square dots. 
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FIG. 2. The scaled exponent B{s) as a function of the effec- 
tive impurity e. The numerical result is shown by the square 
dots. Dashed and solid lines c orres ponds to the analytical 
results for the weak impurities ( |3.19| ) and for the strong im- 
purities (3.28), respectively. 



Analytical expressions of i?(e) can be obtained in the 
limiting cases [Q. For the weak impurity (e <C 1), B{e) 
is evaluated by using the bounce solution in the homoge- 
neous case e = as 



So-e 
31.00 



drVB(0,T';e = 0) + 0(£2) 
16.43e-HO(e2). (3.19) 



This result is shown in Fig. g by the dashed line. 

For the strong impurity £ ^ £ci, the local modulation 
of the field becomes relevant, and a one-mode approxima- 
tion is effective. We first evaluate the solutions of ( |3.S| ) 
for £ = id- The bounce solution <pb{x') agrees with the 
stable solution ipolx'), and is obtained analytically as 



ipB{x';ec\) 



cosh^((2;'-|-aci)/2)' 



(3.20) 



where tanh(aci/2) = 1/^/3- At the critical point e = id, 
there appears the zero frequency mode {X\ — 0) around 
the solution (p-B{x;eci)- The eigenfunction of this mode 
is obtained from (3.12) as 



4''\x';e)^ 



Csinh((a;' -Haci)/2) 
cosh3((x'-t-aci)/2) 



(3.21) 



where C is a normalization factor determined from 
/ \ipi^''\'^dx' — 1, and is given by 



135 



8(9 - 2\/3) 



1.746. 



(3.22) 



When e is slightly smaller than id, this mode describes 
slow dynamics of the system, and is related to the tun- 
neling process. Hence, we can approximate the bounce 
solution as 



(pB{x',T';i) 



(pBix';id) 



X{t')^^^\x'), (3.23) 



where X{t') is a tunneling variable denoting the dynam- 
ics of the local deformation at th e imp urity site. By 



substituting ( 3.23 ) into the action ( 3.17 ), we obtain the 
action in terms of X{t') as 



S[Xir')] = 







(dX\ 


/ dr' 


\ 




'o 




Kdr'J 



V{X) (3.24) 



V{X) 



^iid-i)CX-^iCXf. (3.25) 



This action corresponds to one-particle dynamics in a 
cubic potential V{X). The potential barrier AU and 
frequency of small oscillation around the metastable state 
are given as 



AU = 2^/2 X 3-5/4(£,i - e)3/2 
= 2^/4 ^ ^-ii/sc(^ed - sY'^ 

By applying the standard bounce technique 
this potential, the function B{i) is obtained 



36A[/ 



11.25(eci - e) 



(3.26) 
(3.27) 

M to 



(3.28) 



This result, shown in Fig. H by a solid line, agrees with 
the numerical result for the strong impurity strength, and 
gives a goo d est imate even for the weak impurity. The 
expression ( |3.28 ) may be useful to analyze experimental 
data. 

In order to analyze experimental results, it is conve- 
nient to draw the graph of B{i). The feature of impurity 
effects can be obtained by this analysis. Concerning the 
highly-biased region, this plot is universal in the sense 
that it does not depend on the bare impurity strength e. 

In order to apply the viewpoint of the dimensional 
crossover, the bias-dependence of the exponent B is dis- 
cussed . F or simplicity, the cosine potential with a bias 
term {'Li) is considered with a notation fc,Q = 1, (see 
above Eq. (U)) a = V2, 6 = 1 (see below Eq. <^^). We 
use the bias parameter ry = (/c,o— /)/(/c.o— /c)- It takes 
and 1 corresponding the classical threshold (/ = fcfi = 1) 
in the absence of the impurities and to the critical bias 
(/ = /c) in the presence of the impurities, respectively. 
Here, /c is defined by id = 4/^/3 = £(2(1 -/c))-3/4. The 
graph of i? as a function of 77 is shown in Fig. ^. The 
numerical result is shown by square dots. Analytical re- 
sults for the weak and strong impurities are also shown 
by the dashed and solid curves, respectively. 
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FIG. 3. The scaled exponent B versus the bias parame- 
ter 7] = (/c,o ~ /)/(/c,o ~ /c). The numerical results are 
shown by square dots. Dashed and solid curv es co rrespond 
to the analytical res ult fo r the weak impurity (3.19) and for 
the strong impurity (3.2J), respectively. The horizontal dot- 
ted line shows the homogeneous case e = 0. The solid and 
dashed curves in the inset show the behavior o f B near the 
threshold bias r; = 1, and the asymptotic form (3.29) around 
ri = 1, respectively. 



The dimensional crossover discussed in Sec. [II B can 
be clarified in this figure. In the region r/ ^ 1 where the 
effective impurity strength s becomes small, B becomes 
almost independent of the bias as seen in Fig. H and the 
exponent B behaves as -B cx (fc.o — f)^ as expected in 
(3.14). As 77 approaches the critical threshold rj — 1, the 
effective impurity strength is enhanced, and the exponent 
is suppressed. The analytical expression near 77 = 1 is 
obtained from ( 3.2S| ) as 



B ~ 22.36(ry- 1 



,5/4 



(3.29) 



Hence, the expo nent behaves as B oc (/c — /)^^^ as ex- 
pected in ( 3.15| ). We show the bias-dependences near 
= 1 in the inset of Fig. p|. The exponent, thus, behaves 
as i? oc (/c — /)3/2-('i+iT/4^ g^jjj ^jjg spatial dimension 
d seems to be reduced to zero by the strong impurities 
accompanied with the change of the classical threshold 

fc ^ 

Finally, we show the bounce solution ip-B{x,T) in Fig. 
^ in the region x > 0, r > 0. Note that the bounce so- 
lution has a reversal symmetry ip-B{x,T) = (/3b(— a;,r) = 
(Pb{x,—t). The bounce solution is expected to describe 
the feature of the nucleation process; ip{x,T) describes 
the shape of the field variable ip at the imaginary-time t. 
As the effective impurity strength increases, the bound- 
ary condition for the bounce solution (pBix, ±00) — (Pq{x) 
is modified at the impurity site. As seen in Fig. the 
nucleation process occurs only at the small region near 
the impurity site for the strong impurities. This change 
of the nucleation process is responsible for the qualitative 
difference between the weak and strong impurities. 





FIG. 4. Bounce solution for (a) £ = 0, (b) e = 1.24, 
(c) e = 1.98. The critical value at the threshold bias is 
e = 4/\/3 = 2.31. 



D. Prefactor 

The prefactor A is ordinarily a minor factor which is 
hard to be measured experimentally. However, it de- 
scribes the system-size (L) dependence of the nucleation 
rate, which may be observed. Here, we study the prefac- 
tor A analytically in the limiting cases within the bounce 
method, and discuss the L-dependence of the nucleation 
rate. 



The prefactor is formulated in (3.11) with the improve- 

0). In the ab- 



ment (3.13) due to the zero mode (Aj^'' 



sence of the impurities, not only the second but also the 
third lowest modes become the zero modes (A3 — 0), 
which need to be treated by the Fadeev-Popov method 
as 



'aBjO) 
27r 



L/2 



L/2 



(3.30) 



6 



where x'q denotes a position of the bounce center. This 
treatment gives the correct system-size dependence, T oc 
L. 

As for the weak impurity, the frequency of the third 
lowest mode is hfted to a small posit ive frequency. Also 
in this case, the replacement in ( 3.30 ) is needed in order 
to obtain the correct system-size dependence. As a result, 
the tunneling rate is calculated as 



r = rn 



L 



L/2 



L/2 



dx'Q{e 



1) 



(3.31) 



where Fg is the nucleation rate per unit length in the ab- 
sence of the impurities, and the function /(xq) is defined 
by 



dr' ipb{x'q,t' ;£ = 0). 



(3.32) 



From this expression, the tunneling rate can be divided 
into two parts as 



F 



bulk 



imp — Tq I u-J^o 
'-L/2 



L/2 



dx'^{e''^f^<^ - 1). 



(3.33) 
(3.34) 



The bulk part Fbuik is proportional to the system size 
L, while the impurity part Fi,„p is independent of the 
system size for L ^ 1. 

Within the bounce method, the L-dependence is natu- 
rally described as follows. For large L, Fbuik alway dom- 
inates Finip, and no impurity effects appear. As L is re- 
duced, the bulk part is suppressed, and below a crossover 
length Lcrj the impurity part may overcome the bulk 
part. Only in this region, the impurity effect is clearly 
observed experimentally. The ratio Fimp/Fbuik is deter- 
mined by t he scaled length L and the impurity factor 
ere in ( 3.34 ). For L 3> 1, the crossover length Lcr can 
be obtained as a function of ae by taking the ratio as 



/oo 
dTVB(a; = 0,T';e = 0), 
-OO 



(3.36) 



as log Lcr = AS. Since Ai? is estimated by the local 
suppression of the potential barrier A[/ and the typical 
attempt frequency loq as Ai? ~ AU/hwo, the expression 
in (uM) can be reproduced qualitatively. 




FIG. 5. Crossover length L, 



imp /Fbuik 



1, 



Thus, the phase diagram is obtained by ^2:23 d^-^^ ) 



(/c — f)^^^/b. In the region L > Lcr, nucleation 



occurs in the bulk region, while in the region L < Lcr, the 
impurity contribution Pimp becomes dominant. The solid and 
dashe d lines show the result and, the asymptotic expression 
( |3.35| ) valid in the limit ai ^ 1, respectively. 

Beyond the perturbational analysis of e, it is difhcult to 
calculate the prefactor, because the low-frequency modes 
may couple with each other. It would be a future problem 
to co nside r the nucleation rates in this regime. (See also 
Sec. VB.) Only for the strong impurity, the prefactor 
can be calculated by assuming L <g; Lcr and by applying 
standard procedures p^,E9^,p8| to one-variable problem 



Fig. |5|, where the crossover length is shown by the solid 
curve. In the upper region of the phase diagram, the nu- 
cleation occurs dominantly in the bulk region, while in 
the lower region, the impurity part dominates the bulk 
part, and inhomogeneity effects appear clearly. 

The intuitive discussion about the appearance of the 
impurity effect in Sec. || can be now clarified in de- 
tail. The asymptotic form of the crossover length Lcr for 
ere S> 1 is obtained analytic ally b y applying the station- 
ary method to the integral (3.34) as 



A - V60\ 



2tt ' 



(3.37) 



log Lcr ^ 16.43(cre) - 0.211 - i log(cre). (3.35) 

This result is shown in Fig. |^ by the dashed curve. The 
predominant term in ( p.35| ) can be related to the en- 
hancement of the exponent AB defined by 



although the crossover length Lcr cannot be easily ob- 
tained even in this case. 



IV. THERMAL ACTIVATION REGIME 

At high temperatures, nucleation is caused by thermal 
fluctuations. Since the nucleation rate has already been 
studied in this region in ref. , we only summarize the 
feature of the impurity effect. 

In the thermal activation regime, the nucleation rate 
is evaluated by Kramers' formula 



7 



r = Aexp(-At7/fcBr). 
The potential barrier AU is calculated as 

AU = U[^Bix')]~U[ipoix% 



62 



dx' 



L/2 



1 



1 



,2 „ 3 „ ^S{x'Mx') 



(4.1) 



(4.2) 



(4.3) 



Here, (Pq{x') {cpb{x')) is the stable (unstable) stationary 
solution satisfying 



SU 
Sip{x') 



2 



eS{x') = 0. 



(4.4) 



The same result can be derived by applying the bounce 
method to a finite temperature region p9|-^. Within 
the bounce method, the feature at high temperatures is 
explained as follows. At high temperatures, the bounce 
solution becomes independent of r', since the interval in 
the r direction becomes short. It means that the tem- 
poral dimension become irrelevant. Hence, the relevant 
dimension of the bounce is reduced from (d+l) to d. Ex- 
cept for the absence of the temporal dimension, we can 
discuss the dimensional crossover in the same way as Sec. 
[HB: the energy barrier behaves as AU oc (/c — f)^^'^ for 
the weak impurity, while as AU oc (/c — /)'^^^ for the 
strong impurity. 

The barrier height is calculated in a form 



AU 
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AU{e). 



(4.5) 



The impurity effects appear only through the normalized 
potential barrier AC/. We show AU as a function of the 
effective impurity strength e in Fig. ^ by the square dots. 
The critical impurity strength is given by eci = 4/\/3 
which is the same as the quantum tunneling regime. The 
exponent is obtained analytically for the weak impurity 
e <C 1 as 



AU{e) = 24/5 - Se, 
while for the strong impurity e ~ e^ci as 
AU{e) = 1.433(£ci-e)^/^- 



(4.6) 



(4.7) 



These results are shown in Fig. ||by the dashed and solid 
curves. The feature of the impurity effect is the same as 
the zero-temperature case shown in Fig. 0. 



AU 



T \ \ r 



0. 5 



1. 5 



2. 5 



FIG. 6. The normalized barrier height A[7 versus the effec- 
tive impurity strength e. The numerical results are shown by 
the square dots, while analytical results in the limiting case, 
(4.6) and ([l.7|), are shown by the dashed and solid curves, 
respectively. 



V. DISCUSSION 

A. Extension to the d> 2 case 

First, let us discuss the two-dimensional systems. The 
stable solution ipQ{x',y') is determined by 

fx'x' + fy'y' = - e6{x')S{y'), (5.1) 

where the last term describes the impurity located at 
{x',y') = 0. For the stationary solution having a ro- 
tational symmetry around the origin {x\y') ~ (0,0), 
this equation is reduced, by using the radius coordinate 
r = ^x'2 + 1/'2, to 



I d ( 1 2 . . 



(5.2) 



Here, g{r) is a modified delta function satisfying g{r) = 
for r > 0, g{r) — oo for r = 0, and 



/•oc 

/ dr2TTrg(r) — 1. 
Jo 



(5.3) 



From ( p.2[ ) and ( x2 ) , the boundary condition at r = is 
obtained as 



dip 
iim r— — 

1 — »o or 



s 

2^' 



(5.4) 



Since the function ripr is a continuous function at r > 0, 
if we take an arbitrary constant A (< e/27r), there exists 
a constant a, and for all r in the range < r < a 
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or 



(5.5) 



is satisfied. As a result, ip{r) is evaluated for < r < a 
as 



(^(r) > ^log(^) +</7(a), 



(5.6) 



and the value of Lp diverges at r = 0. Since must be 
finite in the present model, this divergence means that 
the local inhomogeneity cannot be described by the delta 
function in the two-dimensional systems. In other words, 
if we approximate the delta function as g{r) = l/Trrp for 
< r < To and as g{r) = for tq < r, the value of (p(0) 
depends crucially on the impurity size tq, and diverges 
in the limit rg 0. Also in the (i = 3 case, the field 
</?(r) diverges at r = 0, since the bounce solution behaves 
^p{r) (xl/r near r = 0, where r — x'"^ + y''^ + 0'^, 

These results indicate that in the highly-biased region, 
the system is sensitive to local inhomogeneities. To clar- 
ify the inhomogeneity effects in the d > 2 case, it is 
needed to consider the detailed profiles of the inhomo- 
geneities beyond the local approximation adopted in this 
paper. 



B. Validity of the bounce method 

In this paper, the nucleation rate is calculated by the 
bounce method. Although this method has been applied 
to multi-variable tunneling problems for many years, we 
should be careful of its validity. In principle, the exact 
tunneling rate should be determined by the full informa- 
tion about the action S[(j){x')]. The bounce solution can 
extract important informations from the action; the ex- 
ponent i? is a difference of the action between the bounce 
and stable solutions, while the prefactor A is determined 
by the second-order fluctuations around both the solu- 
tions. While the bounce method is a convenient method 
to evaluate the tunneling rate, some informations which 
may affect the tunneling rate are dropped. 

It is expected that the bounce method can be justi- 
fied in the following conditions |38|: (1) we can define 
the tunneling variable describing slow dynamics among 
the many degrees of freedom, and (2) this tunneling vari- 
able is well separated from the other variables with fast 
dynamics. In these conditions, the system can be re- 
duced to the one-variable tunneling problem where the 
bounce method is justified well. Here, let us consider the 
model adopted in this paper. For the strong impurity, 
it has been shown that the single-mode approximation 
gives a good estimate of the exponent B. In this region, 
the bounce method is expected to give correct results for 
both the exponent and the prefactor for L <^ La, though 
the crossover length Lcr itself cannot be evaluated. For 
the moderate impurity, however, there appears another 
low-frequency mode related to the Goldstone mode. This 



low-frequency mode is expected to couple with the tun- 
neling mode. Hence, the bounce method may give an 
incorrect result especially for the prefactor determined 
only by the information about second-order fluctuations. 
It is a future problem to improve the bounce method in 
this situation. 



VI. SUMMARY 

The local inhomogeneity effects on nucleation pro- 
cesses have been studied within the bounce method at 
high biases. The nucleation rate has been calculated in 
a form F = Aexp{—B/h). It has been found that the 
effective impurity strength e can be controlled not only 
by the bare impurity strength e but also by the exter- 
nal bias /. The exponent B has been calculated as a 
function of the effective impurity strength e without re- 
duction to a single-variable tunneling problem. It has 
been shown that these results are well reproduced for 
the weak impurity by a perturbational analysis and for 
the strong impurity by a one-mode approximation. It has 
been clarified that the results can be explained by the re- 
duction of the effective dimensionality of the system due 
to the impurity. By calculating the prefactor A, the con- 
dition of the appearance of the impurity effects has been 
discussed, and the phase diagram has been obtained. 

There remain future problems. One of them is to study 
the inhomogeneity effect in higher-dimensional systems. 
As discussed in Sec. the local approximation does 
not hold for the d > 2 case. From this result, it is 
conjectured that the sample inhomogeneities may be- 
come essentially important for higher-dimensional sys- 
tems. I think that the study of the inhomogeneity effects 
in high-dimensional systems would give an important as- 
pect about the interplay of quantum tunneling and ran- 
domness. 
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APPENDIX A: RELEVANCE OF THE LOCAL 
IMPURITY MODEL TO EXPERIMENTAL 
SYSTEMS 

In this paper, the inhomogeneity effect is described by 
the potential 



(Al) 



If the function e{x) changes only in a narrow region at 
the impurity site compared with the typical length scale 
(~ a^^/^(/ — /c)~^^^), it can be replaced with the delta 



function S{x) in (2.6). The potential form (|Al|) can be 
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realized in various experimental situations. One typical 
example is a spatial inhomogeneity of the bias. If the 
bias is assumed not to be uniform but to have a spa- 
tial de pen dence as / + s{x), then we have the potential 
form ( |A1[ ). Another example is the inhomogeneity of 
the amplitude of the potential barrier in VolcjA^ If we 



replace, for example, the cosine potential in (2.3) with 
(1 -|-e(a;))(l — cos 0), then the form (Al) is also obtained. 
As for the Josephson junction systems, the coupling to 
the derivative of the field (jfxjx) can also be introduced 
by applying the magnetic field [E5| . If the cou pling is de- 
noted with b{x)(l)x{x), the potential form ( |Al[ ) is derived 
by taking the function e{x) as —bx{x). 

In addition to the above situations, boundaries play a 
role of local inhomogeneities. The potential barrier AU 
is reduced to AU/2 at open edge. Hence, the nucleation 
occurs dominantly at the boundary. The prefactor A of 
the nucleation rate is always independent of the system 
size, while the exponent B is affected by the boundary 
condition. Assuming that the field (p{x') is defined at 
x' > 0, and that the boundary condition is given as 



ipx'{x' = 0) = /i, 



(A2) 



it can be shown easily that this boundary effect corre- 
sponds to the local inhomogeneity by taking e — 2h. 
Thus, the exponent -Bodge is obtained as 

Sedge = ^<jB{e = 2h). (A3) 

By changing the boundary condition at the edge, the 
impurity effects on the exponent can also be observed. 
Actually, this tendency has been observed in experiments 
0. 



[1] P. Hanngi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 

62, 251 (1990). 
[2] I. M. Lifshitz, and Yu. Kagan, Zh. Eksp. Teor. Fiz. 62, 

385 (1972) [Sov. Phys. JETP 35, 206 (1972)]. 
[3] T. Satoh, M. Morishita, M. Ogata, A. Savada, and T. 

Kuroda, Physica B 169, 513 (1991). 
[4] T. Satoh, M. Morishita, M. Ogata and S. Katoh, Phys. 

Rev. Lett. 69, 335 (1992). 
[5] S. Balibar, C. Guthmann, H. Lambare, P. Roche, E. Rol- 

ley, and H.J. Maris, J. Low Temp. Phys. 101, 271 (1995). 
[6] J. P. Ruutu, P. J. Hakonen, J. S. Penttila, A. V. Bobkin, 

J. P. Saramaki, and E. B. Sonin, Phys. Rev. B 77, 2514 

(1996). 

[7] A. O. Caldeira, and K. Furuya, J. Phys. C, 21, 1227 
(1988). 

[8] E. M. Chudnovsky, and L. Gunther, Phys. Rev. B 37 
9455 (1988). 

[9] E. M. Chudnovsky, and J. Tejada, Macroscopic Quantum 
Tunneling of the Magnetic Moment (Cambridge Univer- 
sity Press, Cambridge, 1998). 



[10] S. Takeuchi, T. Suzuki, and H. Koizumi, J. Phys. Soc. 

Jpn. 69, 1727 (2000). 
[11] B. V. Petukhov, H. Koizumi, and T. Suzuki, Philos. Mag. 

A 77, 1041 (1998). 
[12] K. Maki, Phys. Rev. B 18, 1641 (1978). 
[13] K. Hida, U. Eckern, Phys. Rev. B 30, 4096, (1984). 
[14] S. Nakaya, and K. Hida, J. Phys. Soc. Jpn. 55, 3768 

(1986) . 

[15] Ji-Min Duan, Phys. Rev. B 48, 4860 (1993). 

[16] K. Maki, Phys. Lett. A 202, 313 (1995) 

[17] N. Hatakenaka, M. Shiobara, Ken-ichi Matsuda, and S. 

Tanda, Phys. Rev. B 57, 2003 (1998). 
[18] S. V. Zaitsev-Zotov, Phys. Rev. Lett. 71, 605 (1993). 
[19] S. V. Zaitsev-Zotov, Phys. Rev. Lett. 72, 587 (1994). 
[20] S. V. Zaitsev-Zotov, G. Remenyi, and P. Monceau, Phys. 

Rev. B 56, 6388 (1997). 
[21] A. Davidson, B. Dueholm, B. Kryger, and N. F. Peder- 

sen, Phys. Rev. Lett. 55, 2059 (1985). 
[22] A. V. Ustinov, T. Doderer, R. P. Huebener, N. F. Peder- 

sen, B. Mayer, and V. A. Oboznov, Phys. Rev. Lett. 69, 

1815 (1992). 

[23] M. G. Castellano, G. Torrioh, C. Cosmelh, A. Costantini, 
F. Chiarello, P. Carelh, G. Rotoh, M. Cirillo, and R. L. 
Kautz, Phys. Rev. B 54, 15417 (1996). 

[24] T. Kato, and M. Imada, J. Phys. Soc. Jpn. 65, 2963 
(1996). 

[25] A. V. Ustinov, B. A. Malomed, and E. Goldobin, Phys. 

Rev. B 60, 1365 (1999). 
[26] For example, see ref. (tel. 
[27] S. Coleman, Phys. Rev. D 15, 2929 (1977). 
[28] C. G. Callan, and S. Coleman, Phys. Rev. D 16, 1762 

(1977). 

[29] S. Coleman, The Use of Instanton in The Whys of Subnu- 
clear Physics, edited by A. Zichichi (Plenum, New York), 
p. 805. 

[30] H. A. Kramers, Physica (Utrecht) 7, 284 (1940). 

[31] M. Yumoto, H. Fukuyama, H. Matsukawa, and N. Na- 

gaosa, J. Phys. Soc. Jpn. 68, 170 (1999). 
[32] M. Yumoto, H. Fukuyama, H. Matsukawa, and N. Na- 

gaosa, Physica B 284-288, 1667 (2000). 
[33] M. Yumoto, H. Fukuyama, H. Matsukawa, and N. Na- 

gaosa, J. Phys. Soc. Jpn. 69, 2953 (2000). 
[34] T. Kato, J. Phys. Soc. Jpn. 69, 2735 (2000). 
[35] D. W. McLaughlin, and A. C. Scott, Phys. Rev. A 18, 

1652 (1978). 

[36] Y. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 
763 (1989). 

[37] B. I. Ivlev, and V. I. Mel'nikov, Phys. Rev. B 36, 6889 

(1987) . 

[38] U. Weiss, Quantum Dissipative Systems (World Scien- 
tific, Singapore, 1993). 

[39] J. S. Langer, Ann. Phys. 41, 108 (1967). 

[40] J. S. Langer, Ann. Phys. 54, 258 (1969). 

[41] I. Aflleck, Phys. Rev. Lett. 46, 388 (1981). 

[42] H. Grabert and U. Weiss, Phys. Rev. Lett. 53, 1787 
(1984). 

[43] A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Theor. 
Fiz. 86, 719 (1984) [Sov. Phys.-JETP 59, 420 (1984)]. 



10 



